Production of hypernuclei in multifragmentation of nuclear spectator matter 
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In peripheral collisions of relativistic heavy ions highly excited spectators containing A hyperons 
can be produced. Such strange spectator matter may undergo a break-up into many fragments 
(multifragmentation) as it is well established for ordinary nuclear systems. We generalize the sta- 
tistical multifragmentation model, previously successfully used for the description of experimental 
data, for the case of hypernuclear systems. We predict relative yields of hypernuclei and the main 
characteristics of such a break-up. We point at a connection of this phenomenon with a liquid-gas 
phase transition in hypermatter. 
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One promising way to study hypernuclei is to use the 
copious production of hyperons in relativistic heavy ion 
collisions. As seen in experiments in the GeV range A 
hyperons originating from the hot participant region have 
a broad rapidity distribution and a considerable fraction 
can be found in the spectator kinematical region. In pe- 
ripheral collisions these A hyperons may be captured by 
excited spectators. Subsequently these spectators may 
decay producing cold and possibly exotic hypernuclei. It 
has been demonstrated experimentally that the cross sec- 
tion of this reaction may be on the order of a few micro- 
barns [2] ■ A similar value can be obtained with UrQMD 
calculations by assuming a coalescence of the produced A 
hyperons with spectators [3|. Indeed there are extensive 
plans to investigate hypernuclei produced in such heavy 
ion reactions using stable as well as exotic beams f^- 

An important feature of these reactions is a correla- 
tion between the excitation energy transferred to the 
spectator matter and the number of high energy par- 
ticles produced as a result of collisions between nucle- 
ons in the participant region (see e.g. [1]), which are 
also related to the strangeness production. As a conse- 
quence the excited spectators break-up into many frag- 
ments (multifragmentation) 0, B 0] j ^^'^ ^^'^ captured A 
will be finally attached to these fragments. According to 
the present understanding, multifragmentation is a rela- 
tively fast process, with a characteristic time around 100 
fm/c, where, nevertheless, a high degree of equilibration 
(chemical equilibrium) is reached [8|. This is a result 
of the strong interaction between the nucleons, and we 
believe that a small admixture of other particles like hy- 
perons interacting with similar cross sections Q will not 
change the general statistical behaviour. Multifragmen- 
tation can also be associated with the liquid-gas phase 
transition in nuclear systems, and thermodynamical fea- 
tures of this transition were experimentally investigated 
As was demonstrated by numerous comparisons of 
theory and experiment, statistical models allow a suc- 
cessful description of experimental data [1, 0, [HI, fl^ . 
In this letter we demonstrate that multifragmentation 
of hypermatter can be described in a similar way, and 
some properties of hypernuclei can be explored by exam- 



ining relative yields of hyperfragments. The calculations 
also show that the usage of exotic heavy ion beams will 
indeed allow to produce hypernuclei with extreme N/Z 
ratios and thus hypernuclei which are not experimentally 
accessible with any other technique. 

For the theoretical description of this process we mod- 
ify the grand-canonical version of the Statistical Multi- 
fragmentation Model (SMM) (^ll^. Previously, this ver- 
sion, with explicit addition of mass, charge, energy and 
momentum conservations, was successfully used for the 
description of multifragmentation in relativistic heavy- 
ion reactions 0, UH [H L Microcanonical extensions of 
statistical models, like IJ], are essential only for a few 
details of the phase transition, but not for the general de- 
scription of the produced fragments. The model assumes 
that a hot nuclear spectator with total mass (baryon) 
number Aq, charge Zq, number of A hyperons Hq, and 
temperature T expands to a low density freeze-out vol- 
ume, where the system is in chemical equilibrium. The 
statistical ensemble includes all break-up channels com- 
posed of nucleons and excited fragments. The primary 
fragments are formed in the freeze-out volume V. We use 
the excluded volume approximation V = Vq + Vj , where 
Vb ~ Aq/po {pq «0.15 fm^'^ is the normal nuclear den- 
sity), and parametrize the free volume Vf = kVq, with 

0,[li|,lil. 

Nuclear clusters in the freeze-out volume are described 
as follows: light fragments with mass number A < A 
are treated as elementary particles with corresponding 
spin and translational degrees of freedom (" nuclear gas" ) . 
Their binding energies were taken from experimental 
data d, [H, [11] . Fragments with A = A are also treated 
as gas particles with table masses, however, some excita- 
tion energy is allowed = AT"^ /eq (eq ~16 MeV is the 
inverse volume level density parameter Q), that reflects 
the presence of excited states in *He, ^H, and ^He nuclei. 
Fragments with A > A are treated as heated liquid drops. 
In this way one can study the nuclear liquid-gas coexis- 
tence of hypermatter in the freeze-out volume. The inter- 
nal free energies of these fragments are parametrized as 
the sum of the bulk (Ff ), the surface (F^), the symme- 
try (F^^fj), the Coulomb (F^^), and the hyper energy 
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Here, H denotes the number of A's. The first three 
terms are written in the standard liquid-drop form 



F?(T) 



and FJ^Ij ^j-{A-H- 2Zfl{A - H). The model pa- 
rameters Wo = 16 MeV, Po = 18 MeV, = 18 MeV and 
7 = 25 MeV were extracted from nuclear phenomenol- 
ogy and provide a good description of multifragmenta- 
tion data 5, 7, JT,, J^]. The Coulomb interaction of the 
fragments is described within the Wigner-Seitz approxi- 
mation, and F^^ is taken as in ref. [5|. 

The new term is the free hyper-energy F^j^ . We as- 
sume that it does not change with temperature, i.e., it 
is determined solely by the binding energy of the hyper- 
fragments. Presently, only few ten masses of single hy- 
pernuclei (mostly light ones) are experimentally estab- 
lished pisl Il6j , and there is only very limited information 
on double hypernuclei available. However, there are some 
theoretical estimations of their masses based on a descrip- 
tion of the available data. One of them is the Samanta 
formula (l7j which suggests a hyper term 
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E'^y^^H- [-10.68 + 48.7/ (A'/')). 
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The contribution proportional to A ^/"^ is motivated by 
calculations of the A binding energy in a potential well 
. In what later follows we use F^j^ — E^^^^ as one of 
the versions for our calculations. We have also explored 
another hyper term (in the following we call it the liquid 
drop hyper term): 



F^yP = (H/A) ■ (-10.68A + 21.27^2/3) 



(3) 



In this formula the binding energy is proportional to 
the fraction of hyperons in the system (H/A). The 
second part represents the volume contribution reduced 
by the surface term and thus resembles a liquid-drop 
parametrization based on the saturation of the nuclear 
interaction. The linear dependence at small H/A is 
in agreement with theoretical predictions for hyper 
matter. We have found that for single hypernuclei with 
A > 8 the average deviation of the liquid-drop masses 
from the experimental masses tabulated in jig is 1.8 
MeV, while using the Samanta formula this deviation is 
worse and amounts to 2.2 MeV. If not mentioned other- 
wise we use the liquid-drop hyper term in the following. 

At this point we should note that the properties of the 
primary fragments may change in the medium compared 
to the vacuum due to the proximity of other fragments. 
Presently, there are evidences that the symmetry energy 
(20I |2]] . [2^ and surface energy of hot fragments in 
multifragmentation may be modified. However, these 
corrections are not large and influence multifragmenta- 
tion in finite nuclear systems very little. They are more 
important for astrophysical applications f23|. Moreover, 




FIG. 1; SMM calculations of multifragmentation of specta- 
tor sources with j4o=200, Zo^SO and ylo=100, Zo=40, with 
number of A hyperons Hq—O, 1, 2. Top panel - average mul- 
tiplicity of produced primary fragments (notations are as in 
the middle panel), middle panel - excitation energy of the 
sources (caloric curve), bottom panel - chemical potential ^ 
responsible for strangeness, as function of temperature T. 



since the binding energies of hypernuclei are mostly not 
known even in vacuum, this problem is naturally included 
in searching for a reliable mass formula for hypernuclei 
within this approach. 

The break-up channels should be selected according to 
their statistical weight. In the Grand Canonics this leads 
to the following average yields of individual fragments: 
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Here ^azh is the ground-state degeneracy factor of 

1 /2 

species {A,Z,H), At = {2TTh^/mNT) ' is the nucleon 
thermal wavelength, tojv «939 MeV is the average nu- 
cleon mass. In our case H « A, therefore, the mass 
difference between nucleons and A can be disregarded 
in this expression. The chemical potentials /i, v, and ^ 
are responsible for the mass (baryon) number, charge, 
and strangeness conservation in the system. They can 
be found from the balance equations: 
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Within this extended SMM approach we have per- 
formed calculations for excited spectator sources which 
can be produced during peripheral relativistic heavy ion 
collisions. Below, we show results for heavy systems 
with ylo=200 and Zq=8Q, and intermediate systems with 
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FIG. 2: Multifragmentation of sources with Ao=100, Zq—AO, 
Ho=2. Top panel - average number of hyperons in fragments 
at temperatures T—4, 5, and 7 MeV. Three bottom panels - 
relative yields of fragments (per event), separately for frag- 
ments containing H=0, 1, or 2 A, at the same temperatures. 



Ao=100 and Aq=50, varying the charge-to-mass ratio 
(Z/A) from 0.4 to 0.5, and Hq from to 2. The first 
source has the mass and Z/A ratio typical for a heavy 
projectile. We expect the most extensive absorption of 
A hyperons by spectators in the case of heavy projectiles 
and targets. It will be possible to increase the strangeness 
production by selecting smaller impact parameters, lead- 
ing to a higher excitation and a smaller mass of the pro- 
jectile spectators. Hovifever, previous analyses [l^l sug- 
gest that their Z/A ratio will be nearly conserved. 

In Fig. [1] we present the average multiplicity of frag- 
ments produced after the break-up of the excited specta- 
tors, the caloric curve for these break-ups, and the chem- 
ical potential f , versus the temperature of the specta- 
tor matter. At small temperatures (T < 3 MeV), the 
compound nucleus channel {M = 1) dominates. In the 
following it should decay by means of evaporation or fis- 
sion processes [1| . One can consider it as a nuclear liquid 
phase. At temperatures T = 4 — 6 MeV the system pro- 
duces mainly several intermediate mass fragments (IMF, 
with A: < A< Ao/3), however, a heavy residue and light 
particles are also possible. In nuclear multifragmentation 
this region is interpreted as a liquid-gas phase coexis- 
tence [3, [l^l ■ Here the caloric curve signals an increasing 
heat capacity, associated with a rapid rise of the excita- 
tion energy. And at high excitations (T > 8 — 9 MeV) 
the system disintegrates mainly into light particles (a gas 
phase). One can see that the small admixture of A's does 
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FIG. 3: Comparison of SMM calculations with the liquid-drop 
and Samanta descriptions of hyper terms in the mass formula, 
for the same sources as in Fig. 2. Top panel - the strangeness 
chemical potential £, versus temperature T. Middle panel - 
average number of A hyperons in fragments, and bottom panel 
- yields of fragments with two A, at T=4 MeV. 
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FIG. 4: Charge distributions of primary fragments with 
yl=30 produced after break-up of sources with Ao=50, Hq — I, 
Zo=20 (solid and dashed lines), and Zo=25 (dotted and dot- 
dashed lines), at temperature T—5 MeV. H = stands for 
fragments without A, ii" = 1 are fragments with one A. 



not change this general behaviour. However, a drop of ^ 
at T > 4 MeV signals a decreasing average strangeness 
content of the fragments. 

In the top panel of Fig. [5] we explicitly demonstrate 
this decrease of the average number of A's contained in 
the fragments with rising temperature. In three bottom 
panels we show how the mass distributions of fragments 
evolve with the temperature. Around the onset of mul- 
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tifragmentation (T « 4 MeV) there is a "U-shaped" dis- 
tribution with a big residue and few small fragments. 
The residue contains nearly all strangeness of the system, 
while the small fragments are mainly without hyperons. 
Also the share of free A's is not essential, hyperons pre- 
fer to be bound in fragments in all cases. When the 
system breaks into many fragments, the strangeness is 
distributed among all produced IMF. Within the phase 
coexistence region (T « 5 MeV) we have a rather flat 
distributions for the produced big residues, and nearly 
equal probabilities for producing the residues with differ- 
ent number of hyperons. Large fluctuations of number of 
fragments and their size are expected in this case. 

We have found that there is a sensitivity of the frag- 
ment yields in multifragmcntation to the mass formulae 
used for description of binding energy of hypernuclei. In 
Fig. [3] we compare SMM calculations performed with the 
liquid-drop hyper term ^ in free energy of individual 
fragments, and with the Samanta term There is 

a clear difference in the chemical potential and, as 
a result, the yields of hyper- fragments are also differ- 
ent. As one can see from the bottom panels the liquid- 
drop formula predicts more strangeness in IMF's than the 
Samanta formula. The difference in the yields is partic- 
ularly large for small double hyper- fragments. In future, 
this observable may allow to test experimentally different 
mass formulas for hypernuclei in multifragmcntation. 

It is an important feature that the isotope composition 
of the produced hypernuclei can be varied considerably 
in these reactions. In Fig. [5] we show charge distribu- 
tions of fragments with A=30 after the break-up of spec- 
tator sources with different isospin. The corresponding 
Zo/Aq ratios can be easily obtained by selecting different 
projectiles. After multifragmcntation, the hot primary 
fragments have approximately the same average isospin 
content as the sources, whereas the widths of the charge 
distributions depend mainly on the fragment symmetry 
energy . Thus the present work shows that it is indeed 



possible to produce very neutron rich hyper-fragments. 
We expect that the final isospin distributions will be sim- 
ilar to the ones reported in previous multifragmentation 
studies with exotic beams [8]. 

After the break-up, the Coulomb acceleration and the 
secondary de-excitation of primary hot fragments should 
be taken into account [Sj . As it is well known from calcu- 
lations and experiments, the de-excitation of nuclei with 
A < 200 will proceed mainly by emission of nucleons. As 
a consequence the resulting distributions of cold nuclei 
are not very different from the primary ones, they are 
just shifted to lower masses. In this case the reported 
regularities will not change, and the differences in rela- 
tive yields will survive, though some of them may become 
smaller [2l|, . The last stage will be the mesonic or 
non-mesonic decay of hypernuclei which can used for the 
identification of the hypernuclei 0. 

We conclude that multifragmentation reactions offer a 
new possibility for investigating hypernuclei under con- 
ditions essentially different from those accessible in con- 
ventional nuclear structure studies. One of the advan- 
tages is that in this reaction one can abundantly produce 
hypernuclei with unusual Z/A ratios. It is encouraging 
that one can distinguish between different mass formu- 
lae of hypernuclei, since their properties are manifested 
in their relative yields. More generally, we suggest that 
this reaction can be considered as an experimental tool to 
study clusterization of nuclear matter with strangeness, 
and the phase diagram of hyper-nuclear matter at den- 
sities p « 0.1 — 0.3po a-nd temperatures around T « 3-8 
MeV reached in the freeze-out volume. Apparently, such 
a phase transition may take place in some astrophysical 
cases, for example, in neutron star crusts. 
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